LOW BANDWIDTH ROBUST CONTROLLERS FOR FLIGHT 

NASA Grant NCC 2-711 


June 1991 - December 1992 


Final Report 


by 


Dr. Daniel J. Biezad 
Principal Investigator 
(805) 756-5126 

and 

Hwei-Lan Chou 
Graduate Student 

Aeronautical Engineering Department 
California Polytechnic State University 
San Luis Obispo, CA 93407 


presented to 


Mr. Glenn Gilyard 
NASA Technical Officer 
(805) 258-3724 


May 1993 


Propulsion and Performance Branch, D-2112 
XRP PO 273, NASA Dryden 
Edwards, CA 93523 



Table of Contents 


Item E a££ 

List of Figures 

Ust of Tables ™ 

Preface „ 

Acknowledgements v ^ > 

Synopsis y** 1 

Abstract “ 

1. Introduction * 

2. Strategy of Throttle-Only Flight Control 3 

2.1 Pitch Control 

2.2 Yaw-Roll Control 

2.3 Speed Control 

2.4 Couple Throttles Command to Stick/ 

Thumbwheel Motion 

3. Boeing 720 Linear Model 4 

4. Engines and Bare Airframe System Survey 4 

4.1 Engines 

42 Bare Airframe 

5. Overview of Quantitative Feedback Theory (QFT) 6 

6. Quantitative Feedback Theory (QFT) Control Design 8 

6.1 System Modeling 

6.2 Performance Specification 

6.3 Airplane Parameter Uncertainty 

6.4 Controller and Prefilter Design 

7. Results and Discussion 14 

8. Conclusions and Recommendations 17 

References 1® 

Appendix A: B-720 Configurations 

Appendix B: Papers Produced in Support of this Grant 



List of Figures 


Figure 1. 
Figure 2. 
Figure 3. 

Figure 4. 

Figure 5. 

Figure 6. 
Figure 7. 
Figure 8. 
Figure 9. 
Figure 10. 
Figure 11. 
Figure 12. 
Figure 13. 
Figure 14. 
Figure 15. 
Figure 16. 
Figure 17. 

Figure 18. 
Figure 19. 


Boeing-720 

Boeing-720 simulation cockpit 

Comparison of the longitudinal open-loop response of the B-720 linear 
and nonlinear model, 20% step throttle cmd, nominal configuration 

Comparison of the lateral open-loop response of the B-720 linear and 
nonlinear model, 5% differential throttle, nominal configuration 

Comparison of the longitudinal open-loop response of the B-720 modified 
linear and nonlinear model, nominal configuration 

Engine spooling block diagrams and Engine 6 to *TC 

Pitch rate to thrust bode 

QFT unit feedback control structure 

Equivalent MISO systems of a 3x3 MIMO system 

Typical close loop tracking specification 

Typical disturbance rejection specification 

U contour construction 

Performance bounds, U contour, and the optimal L 0 
A prefilter is needed to fulfill the performance specification 
Flight-path-angle control block diagram 
Bank-angle control block diagram 
Flight-path-angle control with inner q-loop closed 

Bank-angle control with inner p-loop closed 

Performance specification in time domain for approach and landing of 

B-720 TOFC 


in 



Figure 20. 
Figure 21. 

Figure 22. 
Figure 23. 
Figure 24. 
Figure 25. 

Figure 26. 
Figure 27. 
Figure 28. 
Figure 29. 

Figure 30. 

Figure 31. 

Figure 32. 

Figure 33. 

Figure 34. 

Figure 35. 


Performance specification in frequency domain for approach and landing 
of B-720 TOFC 

Transfer function , its performance bounds B(jw), and U contour on 
Nichols Chart; flight-path-angle control 

Open-loop transfer function, L* , on Nichols Chart 

Frequency plot of the close-loop transfer function with no prefilter 

Frequency plot of the close-loop transfer function with prefilter 

Transfer function G^ , its performance bounds B(jw), and U contour on 
Nichols Chart; bank-angle control 

Open-loop transfer function, L^, on Nichols Chart 

Frequency plot of the close-loop transfer function with no prefilter 

Frequency plot of the close-loop transfer function TJ with prefilter 

B-720 augmented control, step flight-path-angle response with no 
turbulence, nominal configuration 

B-720 augmented control, step bank-angle response with no turbulence, 
nominal configuration 

B-720 augmented control, step flight-path-angle response with 
turbulance, nominal configuration 

B-720 augmented control, step bank-angle response with turbulance, 
nominal configuration 

Flight -path-angle and throttle response to full-forward stick deflection, 
nominal configuration 

Bank-angle and throttle response to full-right stick deflection - QFT 
compensation, nominal configuration 

Response to full-right stick deflection - Heuristic compensation, nominal 
configuration 


iv 



Figure 36. Bank-angle response under turbulance due to c.g. 0 

Figure 37. Bank-angle response under turbulence due to nose boom 0 

Figure 38. Robustness of the flight-path-angle control 
Figure 39. Robustness of the bank-angle control 


V 



List of Tables 


Table 1 . Investigation of Longitudinal Feedback Parameter 
Table 2. Investigation of Lateral Feedback Parameter 

Table 3. QFT Performance Specification 

Table 4. Flight Configurations for B-720 Approach and Landing 
Table 5. Longitudinal Mode Comparison 

Table 6. Lateral Mode Comparison 


vi 



PREFACE 


The statements and conclusions of this report are those of the authors and 
not necessarily those of the California Polytechnic State University, San Luis 
Obispo. The evaluations described in this document are based solely on the 
results of tests conducted by the authors during this project. 

This report does not constitute a standard, regulation or specification. The 
mention of commerical products, their source, or their use in connection with 
materials reported herein is not to be construed as either an actual or implied 
endorsement of such products. 


ACKNOWLEDGEMENTS 

The authors gratefully acknowledge the sponsored research by NASA, 
and in particular acknowledge and thank Mr. Glenn Gilyard and Mr. Bill 
Burcham for their technical support, guidance, and encouragement throughout 
this project. 

The material in this report was prepared and written by Professor Daniel 
Biezad and Graduate Student Hwei-Lan Chou; and edited by ARDFA Program 
Assistant Sue Dotson. 


VII 



SYNOPSIS 


This final report for NASA Grant NCC 2-711 cover reporting period June 
1992 through December 1992. The report analyzes the longitudinal and lateral 
flying qualities of Propulsive-Only Flight Control (POFC) for a Boeing 720 
aircraft model. Using Quantitative Feedback Theory (QFT), performance results 
from compensators are documented and analyzed. This report is also the first 
draft of a graduate thesis to be presented by Hwei-Lan Chou. The final thesis 
document will be presented to NASA when completed later this year. 

The latest landing metrics, related to bandwidth criteria and based on the 
Neal-Smith approach to flying qualities prediction, were used in developing the 
performance criteria for the controllers. The compensator designs were tested on 
the NASA simulator and exhibited adequate performance for piloted flight. 

There was no significant impact of QFT on the performance of POFC in either the 
longitudinal or lateral modes of flight. This was attributed to the physical limits 
of thrust available and the engine rate of response, both of which severely limited 
the available bandwidth of the closed-loop system. 



Abstract 


Through throttle manipulations, engine thrust can be used for emergency 
flight control for multi-engine aircraft. Previous study by NASA Dryden has 
shown the use of throttles for emergency flight control to be very difficult In 
general, manual fly-by-throttle is extremely difficult - with landing almost 
impossible, but control augmentation makes runway landings feasible. Flight 
path control using throttles-only to achieve safe emergency landing for a large jet 
transport airplane, Boeing 720, was investigated using Quantitative Feedback 
Theory (QFT). Results were compared to an augmented control developed in a 
previous simulation study. The control augmentation corrected the 
unsatisfactory open-loop characteristics by increasing system bandwidth and 
damping, but increasing the control bandwidth substantially proved very 
difficult. The augmented pitch control is robust under no or moderate 
turbulence. The augmented roll control is sensitive to configuration changes. 
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1. Introduction 


Through throttle manipulations, engine thrust was found useful in providing 
some controllability for multi-engine aircraft in emergency situations with severe or 
complete flight control system failures (such as hydraulic system failures). Aircraft 
flight control systems are extremely reliable. Current generation aircraft utilize multiple 
and independent control surfaces, hydraulics, sensors, control computers, and control 
cables to achieve a high level of control system redundancy and reliability. Although 
rare, severe flight control system failures do occur. 

NASA Diyden has studied the use of throttles for emergency flight control for a 
range of airplanes 1 *^. Many multi-engine airplanes exhibited some degree of useful 
control capability with the throttles. In general, flying an aircraft in manual mode using 
throttles-only requires a tremendous pilot workload and landing is considered 
extremely difficult to almost impossible. Control augmentation, using feedbacks and 
direct coupling of the throttle command to stick /thumbwheel motion, has greatly 
improved flying qualities, and ground simulation landings can be achieved. 

The primary aim of this current study on Throttles-Only Flight Control (TOFC) is 
to develop an augmented flight path control using throttles-only to achieve safe 
emergency landings. Application of TOFC on a large four-engine jet transport airplane, 
Boeing 720 (B-720) (Figure 1), is investigated. An augmented B-720 TOFC, developed 
and implemented on a high fidelity B-720 flight simulator (Figure 2) by NASA Dryden 2 , 
had obtained good pilot rating by increasing the control bandwidth and the phugoid 
and Dutch-roll damping 2 . 

This report presents an alternative control design technique based on 
Quantitative Feedback Theory (QFT) to further improve the Dutch-roll damping and to 
increase the control bandwidth for better handling qualities. The control design uses a 
linearized B-720 model derived from perturbations of the full non-linear equations of 
motion about trim at an approach and landing flight condition. 

A robust controller is highly desirable for systems with plant parameter 
uncertainty (such as an aircraft undergoing configuration changes). The QFT 
technique 6 - 9 was chosen because it allowed designers to specify a desired close-loop 
response and a performance specification, and then built a controller to meet the 
specification. Most of all, the technique can incorporate plant parameter uncertainty 
and plant disturbances into the control system design by converting them into design 
constraints and then design a controller to have the system satisfy the imposed 
constraints. The controller thus designed guarantees robust performance over full 
range of the plant uncertainty while keeping the disturbance effect to the system 
minimum. 

The desired performance specification may not always be achieved within the 
given control actuation and rate limits. However, the transparency of the QFT 
technique throughout the design process preserves many of the insights which are lost 
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in several of the modem control techniques and thus provides control designers with 
valuable information about the system under investigation. QFT also provides a 
quantitative relationship between the amount of uncertainty and feedback (i.e. the 
magnitude of feedback is determined in proportion to the amount of uncertainty, 
therefore, reduces the possibility of overdesign). 

In this report, the strategy of flight control using throttles-only is introduced. 
The fidelity of the linear B-720 model is examined. An overview of QFT with step-by- 
step procedures is provided, and its application on the design of an augmented flight 
path control using throttles-only for approach and landing of B-720 is presented in a 
summary fashion. Control design results using QFT are compared to the augmented 
control developed in a previous simulation study. 


Nomenclature, Abbreviations, and Acronyms 

(a) short form of (s+a) 

c.g. center of gravity 

C, ^ non dimensional yaw-roll coupling derivative 

C mu non dimensional velocity-pitch coupling derivative 

D. R. Dutch-roll 

G“ ul transfer functions 

K q pitch rate feedback gain 

K r flight path angle feedback gain 

K p sideslip angle feedback gain 

bank angle feedback gain 
q pitch rate (deg/sec) 

QFT Quantitative Feedback Theory 

s.p. short period 

TOFC Throttle-Only Flight Control 

z thrust (lbs) 

6 engine rpm 

6jc stick input (full deflection=l unit) 

y flight path angle (deg) 

0 pitch angle (deg) 

P angle of sideslip (deg) 

<p bank angle (deg) 

c o . natural frequency 

£ damping ratio 

K,wJ short form for S 2 + 2£<y n S + (o D 2 
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2. Strategy of Throttles-Only Flight Control 

The propulsion system of a multi-engine aircraft can be used for heading and 
flight path control. Differential throttles is applied to control roll through yaw, and 
symmetric throttles is applied to control pitch. Speed control by throttles becomes 
ineffective when control systems fail. Other means may be used to change the airplane 
speed as described in below. Throttles are coupled to stick /thumbwheel for easier and 
more conventional control handling. 

2.1 Pitch Control 

Symmetric throttles induces a phugoid mode and a speed change, which in turn 
generates a pitching moment change through speed stability effects-C mu . This is the 
primary source of pitch control. Pitch control may also be generated by other factors 
such as pitching moment change due to thrust line offset, flight path angle change due 
to the vertical component of thrust, and an instant pitching moment change generated 
by engines mounted at different vertical levels, as in the case of B-720. 

2.2 Yaw-Roll Control 

Differential thrust generates sideslip, which in turn generates rolling moment 
changes through wing dihedral and sweep effect-C^. Roll is controlled by applying 
differential throttles to achieve the desired bank angle and thus to make turns and 
heading changes. 

2.3 Speed Control 

Retrimming speed by the use of throttles becomes ineffective when primary 
control surfaces are locked due to control systems failure. When control system failure 
occurs at speeds other than landing speed, retrimming to an acceptable landing speed 
may be accomplished by using other techniques such as lowering flaps (assuming the 
electrically controlled flaps are operative), extending landing gears, moving eg. aft, 
varying stabilizer deflection, or varying the speed between the low and high mounted 
engines. 

2.4 Couple Throttle Command to Stick/Thumbwheel Motion 

Direct coupling of the throttle command to stick/thumbwheel motion has eased 
the pilot's handling of control. The airplane can be controlled in a conventional fashion, 
such as pitch up with stick forward or pitch down with stick aft. 
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3. B-720 Linear Model 

The B-720 linear model is derived from perturbations of the full nonlinear equation 
of motion about trim and is completely decoupled in longitudinal and lateral dynamics. 
All control states and inputs are perturbed independently at the steady state of a 
desired trim condition. The inputs are thrust from each engine. There are four 
configurations given for the study of approach and landing of B-720 TOFC. The state- 
space representation of the linear model of these configurations are listed in Appendix 
A. Of the four configurations shown, configuration 1 is the nominal configuration for 
baseline design. 

The fidelity of the linear model was examined by comparing the open-loop 
response of the linear model with the nonlinear model as shown in Figures 3 and 4. The 
longitudinal response of the linear model was about 30% less in magnitude than the 
nonlinear model, and was, therefore, modified by a correction factor of 1.3 (Figure 3). 
The longitudinal linear model after modification closely portraits the longitudinal 
nonlinear model (Figure 5). Linear design analysis utilized a computer control package 
"Program CC" to assist the design 10 . 

The lateral response depicted in Figure 4 shows that the linear model would closely 
follow the nonlinear model as long as the small perturbation assumption is not violated, 
i.e. the command input should be of a small magnitude and a short duration. Figures 3 
and 4 also illustrate that for the nonlinear model, a flight path angle command would 
induce little coupling in roll /yaw, while a bank angle command would induce 
pronounced pitch coupling. Coupling between longitudinal and lateral modes is 
completely absent for the linear model. 


4. Engines and Bare Airframe System Analysis 

4.1 Engines 

Spool-up and spool-down engine dynamics for the B-720 engine are shown in 
Figure 6a 12 . The empirical transfer function developed is given in short form notation 
by 



275 

(.55X5) 


This equation is illustrated in Figure 6b over low frequency ranges up to 1.0 
rad /sec. 
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4.2 Bare Airframe 


It is apparent from the engine Bode diagram in Figure 6b illustrates that severe 
bandwidth attenuation would occur beyond frequencies of 1 rad /sec. Therefore, it may 
not be possible to increase the closed-loop bandwidth beyond 1 rad /sec within the 
range of available thrust This can be seen in the pitch rate "q" to thrust z transfer 
function, G^ (dee/KCj (refer to Appendix A), of the bare airframe shown in Figure 7. The 

full-order trmisfer function G^ /,<e) shows that 80 dB of gain must be added to yield a 
crossover frequency beyond lTad /sec. This corresponds to 10,000 lbs of full thrust 
from each engine, which is not practical for approach and landing. 


A low order fit to G* d ^ / " e) is also depicted in Figure 7 and is very accurate near 

the phugoid frequency. Piloted flight of the unaugmented aircraft was consistently a 
level 3 2 . The main difficulties were the lightly damped phugoid and the low 
bandwidth throttle control. 

The accuracy of the low order fit near the phugoid frequency means that, to a 
first order approximation 12 , the phugoid frequency and damping are found from the 
following equation: 


2£(o 0 - -X u + 


M u (X n - g) 


a>„ 


M « 


Conventional transport aircraft can be shown to be roughly proportional to M u . 

It should be strongly noted here for the classic case of M u =0 and for negative 
values of M u (Mach tuck) that the aircraft cannot be practically flown with throttle alone 
unless rotational control in pitch is added and difficulties will also be encountered as 
M a becomes small (aft c.g. location). Both of these cases require the addition of an 
effective rotational controller about the pitch axis. This maybe achieved by using 
differential inboard and outboard thrust, provided the inboard engines are a different 
distance from the aircraft xy-plane than the outboard engines. These configuration 
characteristics determine the innate capability for throttles-only piloted control. 


5 



5. Overview of Quantitative Feedback Theoiy 

QFT is a frequency domain control technique that uses a fairly straightforward 
and transparent design approach**. To apply QFT, systems are usually modeled in a 
unit feedback form (Figure 8) where all blocks may present scalar (SISO) or matrix 
(MIMO) system transfer functions. For MIMO systems, a m x m MIMO system can be 
converted into am 2 - equivalent multi-input single-output (MISO) loops (Figure 9). 

QFT techniques allow designers to specify a desired performance specification with 
performance tolerance and then incorporate the tolerance with the plant uncertainty 
and system disturbances to form the design constraints: the performance bounds and 
the U contour (Figures 10, 11 and 12). 

The design constraints are then placed on a Nichols Chart together with the 
nominal plant transfer function, Po- A controller will be selected to reshape P Q to form 
Lo (the nominal open-loop transfer function) as to have Lo satisfying all the design 
constraints of performance bounds and U contour (Figure 13). By having Lo satisfy all 
the design constraints, if possible within the given control actuation and rate limits, the 
system is guaranteed robust over the full range of plant uncertainty. However, the 
system may not completely meet the performance specification (Figure 14). A prefilter 
is usually required to further reshape the system to fully meet the specification. The 
prefilter design is implemented on a Bode plot. 

The basic design procedures of the QFT technique for minimum phase systems 
are accomplished by the following four steps: 

1) Model the system in a unit feedback form to apply QFT. A m x m MIMO system can 
be converted into a m 2 -equivalent MISO system and the coupling between loops can 
be considered disturbance input (Figures 8 and 9). 

2) Specify the desired close-loop frequency response performance specification. Figure 
10 shows the construction of a desired close-loop performance specification with an 
upper bound, Bu ; a lower bound, Bl; a tolerance, d„,; and a maximum peak 
magnitude, M m . 

The tolerance, <5„, is specified to obtain robust performance, and the maximum 
M m is specified to obtain a desired system damping. The upper bound is generally 
synthesized by an underdamped second order dose-loop transfer function (T.F.), 
T u (s) and the lower bound by an overdamped dose-loop T.F., T L (s) with figures of 
merit such as settling time, rise time, peak overshoot or damping ratio, and natural 
frequency, etc. 

A desired disturbance performance specification (Figure 11) needs only an upper 
bound to confine the disturbances. The objective of the technique is to design a 
controller such that the variation of the response due to plant uncertainty lies within 
the specified boundaries and the effect of disturbance is minimized, that is to have: 
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Bl(o))s|TrG 0>)j s Bu(co) 

[TdiG to)[ s |MdQ o>)j 
[TdzG «>J*|MdG «)| 

3) Convert the performance tolerance, <5„, and the maximum Mm, onto Nichols Chart to 
form the design constraints: the performance bounds, Bo(ja>i), and the U contour. 

i) Performance bounds are curves on the Nichols Chart that are determined by 
matching the magnitude of the range of plant uncertainty with the magnitude of 
the performance tolerance, <5„. Therefore, satisfying this constraint guarantees 
the variation of the system response due to plant uncertainties will be no greater 
than <5„. There is a performance bound for each frequency. 

ii) On the Nichols Chart, the U contour is a M-cirde that has the magnitude of M m , 
with part of the circle stretched for uncertainty at high frequencies (same as the 
length V shown in Figure 12). By having the open-loop response not penetrating 
the U contour, the system's damping will be guaranteed no less than the 
damping correlating to M m . The construction of a U contour is shown in Figure 

12 . 

4) Reshape the nominal plant transfer function, Po. Gain /pole /zero compensation 

may be placed on Po to reshape it to satisfy the design constraints. After reshaping, 
Po becomes Lo, and the compensation chosen forms the controller, G 0 as can be 
depicted from the relationship: 

L 0 =P 0 *G C 

To satisfy the design constraints, Lo should not penetrate the U contour, while each 
frequency u>i on L 0 should be kept on and above its corresponding BofltOi). The U 
contour , the performance bounds and the optimal Lo of an example problem are 
shown in Figure 13. 

5) After reshaping, the system is guaranteed robust over the full range of plant 
uncertainty, i.e. &r(jwi) s 6 R (jWi) (Figure 14). However, the system may not have 
met the performance specifications completely. A prefilter is usually required to 
further reshape the system to fully meet the specification. 
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6. Quantitative Feedback Theory Control Design 

A QFT computer control package was used to assist the QFT design 8 . The 
program is to be used for minimum-phase plants only, i.e., the plants should have no 
zeros in the right half s-plane, therefore, only the gain curve of the desired performance 
will be specified and satisfied. For nonminimum-phase plant, the phase of the desired 
dose-loop performance shall also be specified and satisfied. 

6.1 System Modeling 

The block diagrams of flight-path-angle control and bank-angle control are 
presented in Figures 15 and 16. The inner pitch-rate (q) loop and sideslip-angle (fi) loop 
were first closed with G^‘ =1, K q =60 and G^ =10, K,=4, respectively, which were the 

heuristic settings chosen by investigating the properties of the inner loop. Tables 1 and 
2 summarize the investigation: 


Table 1. Investigation of Longitudinal Feedback Parameter 


Feedback 

Parameter 

Phugoid Mode J 


<On 

q 

Increase 

(Require high gain) 

No change 

Y 

No change 

Increase 


Table 2. Investigation of Lateral Feedback Parameter 


Feedback 

Parameter 

Lateral Phugoid 

Dutch Roll | 




<»n 

P 

Increase 

Increase 

Decrease 


r 

Small 

increase 

No change 

Small 

increase 

No change 

p 

Small 

decrease 

Small 

increase 

Increase 

No change 

<p 

Increase 

Increase 

Small 

increase 

No change 


To apply QFT, with the inner loop closed, the outer loops are rearranged in a unit 
feedback form as shown in Figures 17 and 18. 


































6.2 Performance Specification 

To obtain good handling qualities, the close-loop response for each of the y- and 
<p- loops, which are also the pilot control open-loops, should have the foUowmg 
characteristics: 

1) A bandwidth, aw * 2 rad / sec for landin & of a transport aircraft 11 . 

2) A k/s gain curve slope (-20 dB/decade) around the crossover frequency, co c “. 

A desired close-loop specification was synthesized based upon these two 
requirements, and is shown in Figure 19, which has a k/s slope near to c - 1.5 rad/sec 
(with 0)o>b - 2 rad / sec) and a comer frequency, co CO mer - 0 8 rad / sec. 

The desired close-loop specification is synthesized in the following four steps: 

1) Synthesize the initial B u and B L . B u is usually modeled by an underdamped 
simple second order close-loop T.F., 




+ 2gnr n S + vj„ 


2 ' 


while B l modeled by an overdamped simple second order close-loop T.F., 


J , whereOj and o 2 s n r n . 

1 (S + C>iX S + a 2> 

With a desired performance specification of g =.6 and m n =.8 rad /sec, this 
yields: 


T (S) - ■ ; 64(S+ H- and T l (s) - 

'uv*^ -2 . PLCr . AA LV 7 


.385 


S 2 +.96S+.64 ' (S+.55XS+.7) 

2) Add a pole to T L (S) to widen the 6„ at high frequencies. This yields: 


.77 

(S+.55XS+.7XS + 2) 


This is required by the Bode derived theorem which states that ^SJcIgt - 0 , 

i.e. the reduction in sensitivity Sj at the lower frequencies must be 
compensated by an increase in sensitivity at the higher frequencies. 
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3) Add a zero at 1 rad /sec to increase the gain slope from -40 dB/decade to -20 
dB / decade. This yields: 


T„(s) - 


.64(S + 1) 
S 2 +.96S+.64 


and T l (s)- 


.77(S + 1) 

(S+.55XS+.7XS + 2) 


4) Raise the whole synthesized gain curve until a (Do*, <• 2 rad /sec is obtained. 


The magnitude of B u , BLand <5„ at each frequency can thus be obtained and are 
shown in Table 3. 


Table 3. QFT performance specification 


Frequency(rad / sec) 

0.1 

0.3 

0.5 

0.7 

1.0 

ECU 

5.0 

Bu(dB) 

17.0 

17.0 

17.3 

16.0 

13.0 

2.0 

-13.0 

BL(dB) 

16.8 

15.0 

12.3 

9.9 


-7.1 

-23.0 

<5„ (dB) 

0.2 

2 

5 


8.4 

9.1 

15.0 


6.3 Airplane Parameter Uncertainty 

Four configurations are provided for the study of approach and landing of B-720 
throttles-only flight control. The flight condition of these configurations are 
summarized in Table 4. Configuration 1 is the nominal configuration for baseline 
design. 


Table 4. Flight Configurations for B-720 A; 

pproach and Landing 

Config. 

Number 

Weight 

(lbs) 

Altitude 
(Ft MSL) 

Airspeed 

(Knots) 

Flaps 

(%) 

Gear 

up /down 

1 

140,000 

4,000 

160 

0 

up 

2 

140,000 

4,000 

145 

30 

up 

3 

160,000 

4,000 

175 

0 

up 

4 

140,000 

4,000 

155 

30 

JiE 
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A plant transfer function with parameter uncertainty is usually described in a 
maximum and minimum format in order to form the plant uncertainty template, which 
will then be used to determine the performance bounds constraint An example of a 
plant with parameter uncertainty described in a maximum and minimum format is 
shown below: 


Example 

For a plant transfer function 

K 2 

Gfsi . where the parameter variations are: 1 < k < 10 and 1 < a < 10 

1 (s + a) 


G(S) min.-— ^ and G(S) max. - 

5 + 1 S + 1UU 


For y- and cp- feedback loops, the minimum and maximum values of the transfer 
functions GJ ( and G^ f determined from the four given configurations, are shown 
below: 

For v-feedback loop: 

The G£ (d ^ g) of the nominal configuration(config. 1) is: 

G r(de8) confie 1 = .01(.203)[.37 > 3.01] 

& (.562) [.624, .111] [.441, L57] (5.25) 


and the min. and max. Gj£^) are: 


min -' 


.0053 (.162) [.35 , 3.01] 

(.40) [.42 , L48] [.66 , .0l] (5.19) 
.01 (.28) [.46, 3.43] 


G r(deg) max V*"' LZTI 

max - ( . 58) [ 45 > L57 ] [.92 , .14] (5.24) 
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For Co- feedback loop: 


The Gt (d ^i N of the nominal configuration(config. 1) is: 

Pin 


.09 [.47 , 3.65] 


nommal (.gg) [. 81 , .is] [.26 , L07] (5.02) 
and the min. and max. of are: 


.06 f. 45 , 3.65] 


mm ’ (.98) [.60 , .15] [.24 , .93] (5.01) 


.09 r»61 , 4.33] 


G Ji < (£g) max * ( L 03) [L 0 , . 20] [. 29 , L 09] (5. 02) 


The QFT control package, used to assist the design, allows the designer to input 
plant parameter variations by entering the transfer function's maximum and minimum 
values for gain, first order poles and zeros, and second order poles and zeros. The 
program forms the plant uncertainty template with the given maximum and minimum 
values, then uses its CAD capability to graphically determine the performance bounds 
required for the design. There are tradeoffs between plant parameter uncertainty and 
system performances. The wider the spread of the parameter uncertainty, the more 
restricted the constraints; consequently more compensation is required. Therefore, the 
performance specification may need to be relaxed when there is not enough control 
power to provide all the compensation that is required. 


6.4 Controller and Prefilter Design 

Pole /zero /gain compensation maybe required to reshape the plant transfer 
functions of the y- and <p- feedback loops, G£ w and G^, and to satisfy performance 

bounds and U contour constraints. On a Nichols Chart, adding a gain will raise the 
transfer function curve, while a zero will bend the curve to the right, and a pole will 
bend the curve to the left The compensation chosen forms the controller, G c . After 
reshaping G^ and G^ become, respectively, L 7 ^ and (the open-loop transfer 

functions of the y- and <p-loops), where L 7 er - G®* * G^ and - G^J * G^ . Each 
frequency Wj on L 7 r and should be placed on and above its corresponding 
performance bounds, Bo(jcoi), to assure robust performance. In addition, L r er and 
must not penetrate the U contour in order to obtain the desired damping. 

Longitudinal Flight-Path-Angle Control: Transfer function G^ and its performance 

bounds, Bo(j(Oi), and U contour are displayed on a Nichols Chart in Figure 21. All 
frequency points on G^ are below their corresponding Bo(ja>i), hence reshaping is 
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required (Figure 21). A pure gain compensator, G C =G* 1 " =16, raises the curve to just 
touching the U contour (Figure 22). Several lead compensators were tried to further 
reshape the G£ to satisfy all Bo(ju)i) while not penetrating the U contour. The lead 
compensators tried had increased the bandwidth and robustness; however, they also 

reduces the output y to a very small value(e.g., an output y=2 degree for a full stick 
input). Therefore, only the pure gain of 1 6 is chosen as the compensator, this left the 

Bo(ja>i) unsatisfied. The frequency response of the dose-loop transfer function, 

(where = L r , r / (1 + L r er ) - (G *G'J/(1 + G£ * G^ ), is shown in Figure 23. It can 

be seen in Figure 23 that <5 T (the spread between Tmax Tmin) l^d exceeded the <5„ 
over the frequency range 0.1 to 0.7 rad / sec as a result of L^, not satisfying the 
performance bounds over that frequency range. To have any frequency, (o u on L^ r 
higher than its corresponding Bo(ja)j) will result in d T (j<u,) ^ while lower than 

Bo(jo)i) will result in <5 T (j<y,) * <5 R (j <*>,)• As can be depicted in figure 23, further 
modification is required to fully meet the prescribed specification. A prefilter of a pure 
gain of 6.5 proved most effective in increasing the bandwidth and met the prescribed 
specification. The frequency response after the prefilter is applied is shown in Figure 
24. 

Lateral Bank-Angle Control: Transfer function G^ and its performance bounds, 

Bo(jcoi), and U contour are displayed on a Nichols Chart in Figure 25. Notice that G^ is 
not only below all performance bounds Boflwj) but it also penetrates the U contour. 
Therefore, more than just a pure gain is required to reshape . A controller, 

Gjp =(s+.15)/(s+1.5), was added to G^ to reshape it and prevent it from penetrating the 

U contour, but it was not successful in satisfying all of the Bo(ja)i). After reshaping, 
is shown on a Nichols Chart in Figure 26. The frequency plot of the close-loop transfer 
function, , where =L^/(1 + L^)-(G^ *G' b )/(G£ *G^>, with no prefilter 
applied yet is shown in Figure 27. A lead compensator of (S+l)/(S+2) is added to 
haunch up the severely deteriorated curve at frequency over 1 rad /sec and to increase 
the phase margin. A lag compensator of (S+0.25)/(S+0.15) is added to steepen the gain 
curve at low frequencies and to provide a smoother k/s curve for good pilot handling 
qualities. The close-loop response after adding the prefilter is shown in Figure 28 and 
the prefilter selected is 15(S+0.25)(S+l)/((S+0.15)(S+2)). 
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7. Results and Discussion 


The objective of this study is to improve the handling qualities for the approach 
and landing of B-720 TOFC by increasing the control bandwidth and the light Dutch- 
roll damping. The control bandwidth of TOFC depends primarily on the engine 
response to throttle command, and on the propulsion-induced low-frequency speed and 
dihedral stability effects, which are configuration-dependent, thus are fixed and 
unalterable. Therefore, the control engineer's only tools are compensation and 
feedback. 


To improve control bandwidth is very difficult as can be depicted from the pitch 
rate to thrust bode in Figure 7. Full thrust of 10,000 lbs from each engine is required to 
yield a crossover frequency just beyond 1 rad /sec. This clearly shows how control 
bandwidth is limited by the control power (the engine) available. 

For flight-path-angle control, pitch rate feedback was effective in increasing 
phugoid damping while y feedback was effective in improving frequency of the 
phugoid mode. For bank-angle control, feedback was found most effective in 

increasing Dutch-roll damping while cp feedback is crucial to lateral phugoid damping. 
Yaw rate feedback, which is effective in damping Dutch-roll when rudder power is 
available, helps Dutch-roll damping and the lateral phugoid damping very little. Tables 
5 and 6 compare the dynamic modes of the bare airframe with the dynamic modes of 
previous simulation designs, of QFT design and of heuristic design(heuristic 

compensation will be discussed later on page 19.) Transfer functions of y to stick and <p 
to stick for all the four configurations are listed in Appendix A. 


Table 5. Longitudinal Mode Comparison 

Density I Phugoid I Short Engine I Pre- G* ta I G* Tc K K 0 

Period filter r ' 

Bare (1.4E-6) (.04,.13) (.65,1.4) (.55)(5) - .... 

Airframe 

Simulation ““ 

Augmented (4.7E-6) (.52,.24) (.52,1.5) (.4)(5.2) 10 1 10 1 4 

Control 
QFT 

Augmented (3.4E-6) (.62,.32) (.46,1.6) (.3)(5.2) 6.5 16 1 1 60 

Control 


//. ' 
i ■ vr i 
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Table 6. Lateral Mode Comparison 




Dutch 

Roll 

Roll 

Engine 

Pre-filter 

Gf" 

G"‘ 


■<p 

K, 

K 

Bare 

Airframe 

(1.1E-4) 

(.12,. 99) 

(1) 

(.55)(5) 

- 

- 

- 

- 

- 

- 

- 

Simulation 

Augmented 

Control 

(.73,. 35) 

(.15,-99) 

(1) 

(5) 

40 

1 

1 

.5 

.5 

1 

- 

QFT 

Augmented 

Control 

(.39) 

(.29,1.0) 

(1.5) 

(.45X5) 

2-5025X1) 

(•i? 

1 

1 

- 

4 

- 


(.15X2) 

(1-5) 






Heuristic 

Augmented 

Control 

(.75,.28) 

(.22,1.0) 

(-9) 

(5) 

40 

1 

1 

.15 

.5 

3 

l 


For longitudinal control, pure gain compensation was used. Since the short 
period mode has a frequency around 1.5 rad /sec (which was beyond the frequency that 
throttles can control) the primary concern was to increase phugoid damping and 
frequency. The phugoid damping and frequency increased from 0.52 to 0.62 and from 
0.24 rad /sec to 0.32 rad /sec, respectively. This increase of response frequency can also 
be depicted from the flight path angle response shown in Figure 29. 

For lateral control, pole /zero compensation was used. The Dutch-roll damping 
was almost doubled, from 0.15 to 0.29. The simulation augmented control has a lateral 
phugoid mode [0.73, 0.35] which combines the spiral and the slow engine mode. This 
was replaced using the QFT design with two real root modes, (0.39) and (0.45), both 
with higher frequencies, therefore faster responses. The comparison of the responses is 
shown in Figure 30. 

All plots in Figures 21 through 32 were obtained from nonlinear simulation runs 
at approach and landing conditions with major control surfaces (ailerons, elevator sand 
rudders) locked while the electrical and mechanical systems, and the landing gear 
remained operative. Figure 33 through 35 shows the throttle response, and flight-path- 
angle and bank-angle fracking response to full stick deflection. 

Turbulence Response: The response of the flight-path-angle control under 
intermediate turbulence is presented in Figure 31. Because of gust randomness, more 
than one simulation run was made to examine the tracking integrity under turbulence. 

The bank angle tracking by QFT design does not perform well under 
intermediate turbulence (Figure 32). This could be caused by the larger Kp gain (Kp=4) 
being used in the p -feedback loop by QFT design, while a Kp=l is used for the 
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simulation designed control. This larger Kp multiplies the gust-induced-sideslip four 
times before it was fed back to the airplane. This had a dramatic effect on the bank 
angle output due to 0.2 degree of sideslip angle would generate approximately 10 

degrees of bank angle, owing to the large Qp of B-720. The sideslip angle (p ) feedback 
is the only parameter that can effectively increase Dutch-roll damping for B-720 TOFC 
A compromise seems necessary between lateral bank angle tracking and Dutch-roll 
damping. 

Good Dutch-roll damping is associated with disturbance excitation of the lateral 
phugoid mode and results in poor tracking of bank angle. Lowering the feedback gain 
reduces Dutch-roll damping but also decreases the sensitivity of bank angle to 
disturbances, and thus makes the lateral phugoid mode less troublesome to the pilot 
when flying in turbulence. 

During the investigation, it was found that the fi being fed back into the B-720 
simulator was the 0 at the c.g. instead of the p at the nose boom. The nose boom p is 

actually measured and fed back into a real airplane. The nose boom p was then 
modeled into the B-720 simulator and the results of the bank angle tracking under 
turbulence were fairly good. Figure 36 shows the bank-angle tracking under turbulence 

due to c.g. p, while Figure 37 shows two runs of bank-angle tracking under turbulence 
due to nose boom p. The p at the nose boom has two more terms caused by lateral and 
longitudinal offsets from the c.g. The dominant term is a function of roll rate. When 
this extra term was active in the feedback loop, lateral performance improved. 

However, there was some question concerning the correctness of the sign of the yaw 
rate term as implemented in the simulation The effectiveness of yaw rate feedback for 
TOFC need to be further investigated. 

A augmented control scheme heuristically determined that feedback p,r, p and 

<p is investigated. The compensation of this control scheme is shown in Table 6. The 
yaw rate feedback is included in the control to improve bankangle tracking, the roll rate 
and bank angle feedbacks are included to increase the damping and frequency of the 

lateral phugoid mode while the p feedback to increase the Dutch-roll damping. The 
heuristic augmented control has a slightly slower response speed ( to n =0.28 < tOn=0.35 of 
simulation augmented control), but a higher Dutch-roll damping (£=.22 > £=.147 of 
simulation augmented control) which has successfully damped the Dutch-roll oscillation 
as can be seen in Figure 35. 

The system response to configuration variations for y-control and for <p-control 
are shown in Figures 38 and 39, respectively. The robustness of the flight-path-angle 
control is improved by QFT as shown in Figure 38. The Dutch-roll oscillation in the 
original simulation compensation is taken out by QFT compensation; however, the 
tracking was not improved (Figure 39). Among the three augmented control developed, 
the heuristically determined augmented control presented the best robust performance 
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(with fairly good bank-angle tracking and no Dutch-roll oscillation.) Due to time 
constraints, the heuristic augmented control was not tested on the real-time nonlinear 
B-720 simulator. 


8. Conclusions and Recommendations 

Studies by NASA Dryden has shown that throttles can be used for emergency 
flight control. Manual fly-by-throttle is extremely difficult with landing almost 
impossible, but with control augmentation, runway landing is feasible. 

Flight path control design using throttles-only to achieve safe emergency landing 
for a transport airplane, Boeing 720, was investigated. Augmented throttles-only flight 
path control built in a previous simulation study has made successful simulation 
landings. However, it showed light Dutch-roll damping and low control bandwidth. 

To increase the control bandwidth substantially proved very difficult. Differential 
throttles to engines mounted at different vertical levels to generate an instant pitching 
moment may be an effective way to increase the control bandwidth. 

For throttles-only pitch control using QFT, the control bandwidth , tracking and 
control robustness were improved by QFT. For bank angle control, QFT has improved 
the Dutch-roll oscillation. However, the lateral phugoid becomes sensitive to 
configuration changes. A compromise is required between Dutch-roll and lateral 
phugoid damping given limited control power. Further investigation of the effects of 
yaw rate feedback is recommended. 
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Figure 1. Boeing-720 
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Figure 2. Boeing-720 simulation cockpit 
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Figure 3. Comparison of the longitudinal open-loop response of the B-720 

linear and nonlinear model, 20% step throttle cmd, nominal configuration 
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Figure 4. 
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Comparison of the lateral open-loop response of the B-720 linear model and 
nonlinear model, 5% differential throttle, nominal configuration 
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Figure 5. Comparison of the longitudinal open-loop response of the modified B-720 linear 
and nonlinear model, nominal configuration 
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Figure 7. Pitch rate to thrust bode 










Dt Da 



where 

R ; Command Step input 
D : Disturbance inputs 
C : output 
P : Plwit 
Oc : Compensator 
F : Prefilter 


The open loop transmission function, L, is defined as: 
L«Gc*P and UnL-LmGc + LmP 


The close loop transfer functions are: 


Tracking: 

Tr- 

FL 

1 + L 

(D1-D2-0) 


Disturbance: 

Tdi - 

P 

(R-02-0) 


1 + L 

Disturbance: 

TD2 - 

1 

(R*0t«0) 


1 + L 


Lm T« * Lm F + lm ■■ - 

1+ L 


Figure 8. QFT unit feedback control structure 




Figure 9. Equivalent MISO systems of a 3x3 MIMO system 
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Figure 10. Typical close loop tracking specification 
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Figure 1 1 . Typical disturbance rejection specification 
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Figure 12. U contour construction 
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Figure 13. Performance bounds, U contour, and optimal Lq 
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Figure 1 5. Flight-path-angle control block diagram 
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Figure 17. Flight -path-angle control with inner q-loop closed 
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Figure 18. Bank-angle control with inner 0-loop closed 



Figure 19. Performance specification in time domain for approach and landing of 

B-720TOFC 
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Figure 20. Performance specification in frequency domain for approach and landing of B- 
720TOFC 
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Figure 22. Open-loop transfer function, , on Nichols Chart 
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Figure 23. Frequency plot of the close-loop transfer function T r r with no prefilter 
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Figure 26. Open-loop transfer function, 
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Figure 28. Frequency plot of the close-loop transfer function with prefilter 
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Figure 29. 
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B-720 augmented control, step flight-path-angle response with no turbulence, 
nominal configuration 
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Lateral Response Simulation compensation 
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Figure 30. B-720 augmentedcontrol, step bank-angle response with no turbulence, nominal 

configuration 
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b) QFT Compensation 

Figure 31 . B-720 augmented control, step flight-path-angle response with turbulance, 

nominal configuration 
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Figure 32. B-720 augmented control, step bank-angle response with turbulance, nominal 

configuration 
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Figure 35. Response to full-right stick deflection - Heuristic compensation, nominal 
configuration 


47 



















Angle of 
sideslip at 
c.g. (deg) 


Sideslip 

angle 

from 

gust(deg) 


TIME TIME.2 
6.000 -[ ~ 

4.000 " 

2.000 - ;°V ’ ‘ 

0.000 “ ^ 

- 2.000 - 

-4.000 - r 

o 0 20.0 

TIME TIME.2 



W*' 


-4.000 - , i 

0.0 20.0 
TIME TiME.2 


Angle of 

6.000 - 

4.000 ” 

2.000 - 

sideslip at 

0.000 ' 

nose boom 
(deg) 

-2.000 - 


-4.000 * 
0.0 



i 20.0 

TIME TIME.2 


igure 37. Bank-angle response under turbulence due to nose boom p 

9 







a) Simulation Compensation 



b) QFT Compensation 


Figure 38. Robustness of the flight-path-angle control 
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b) QFT Compensation 
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c) Heuristic compensation 

Figure 39. Robustness of the bank-angle control 
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APPENDIX A: B-720 CONFIGURATIONS 


The B-720 piloted simulation can be represented by the following block 
diagram: 



Prefilter 


Bank angle loop 
Compensation 


Sideslip Angle Loop ^ 
Compensation . 







The "AIRCRAFT' in the box above represents both the engine and the 
bare airframe dynamics. The engine is approximated by a transfer function, 
Geng, and the bare airframe dynamics are represented mathematically by a 
single quadruple, Pa/c, shown as follows: 


H -n > 



Longitudinal Dynamics 


Throttle command 
% rpm 

< 5 Te (%) 


Engine 


Aircraft Transfer Matrix 



r*i 


a;_b 

throttle inputs equal 

A 

Bi (column)' 

M 

“ 

c ! D J 


C 

D-0 


x -[q( deg/ sec) j a(deg) j v(kts) | 0(deg) j h(ft)] 

y -[Qpiiot * I i q I a I v 1 6 I h I y( de e3 

® - [Zouttxl J lbs ) j Z inbd icftGbs) i z inM right (lbs) | z^ ngh 0bs)] 
Uj - z(lbs) [used wheD all throttles have same command ] 


Lateral Dynamics 


xi rA J B' 

f?l 


A|_ B 2 ' 

x ' 

yj “ [c id 

l u J 


cTd-o 

U 2. 


where 


P_, - 


b 2 

*d = o 


for four engine inputs, u 


for one total engine input, u 2 


x - [p(deg/ sec) j r(deg/sec) j /J(deg) \ ^(deg)] 


y- A y ! p(deg/sec) jr(deg/sec) 1/3 (deg) !^(deg) 

I | | I I J 

® ™ ^outbd left (lbs) j 

^mbd left (lbs) | z^^lbs) J 

^outbd right (lbs)] 

u 2 - [z(lbs)] , where z = z^,* + z inbdkft + (-z^^) + (-z^^) 


The B matrix has four columns, each column is to be multiplied by the thrust 
input from each engine that is given in matrix u. If symmetric throttle is 
given (assume all four throttles are given the same command), the B matrix 
in longitudinal dynamics becomes a single column. Each row value in this 
column matrix Bi is equal to the sum of the corresponding row elements in 




the full order B matrix representing four engines. If differential throttle is 
given (i.e., the left engines and right engines are given same amount of 
command but in opposite directions), the B matrix in lateral djmarrucs 
becomes another single column matrix, B z Each row element mB,is the sum 
of the positive value of columns 1 and 2, and the negative value of columns 3 

and 4 of each row in B. The open-loop configuration then becomes 

PsPa/cTeng, where Peng is the quadruple form of the engine transfer function. 
The quadruples for four different configurations were obtained as 
described in reference. 


Flight 

table. 


conditions for each configuration are summarized in the following 


Configuration Summary 


Config. 

Number 

Weight 

(lbs) 

Altitude 
(Ft MSL) 

Airspeed 

(Knots) 

Flaps 

(%> 

Gear 

up /down 

1 1 

140,000 | 

WKIE22MM 

160 

0 

U P -- 

2 

140,000 

4,000 

145 

30 I 

U -P 

3 1 


4,000 

175 

0 

U P 

4 

140,000 

4,000 

155 

30 

U P 


The transfer functions were obtained from the quadruples using 
System Technology's Program CC. These aircraft transfer functions are listed 
here with each respective row of numbers designating the corresponding 
configuration transfer function values. The nominal configuration, number 
1, is represented by values in each row 1 below. 

l ongitudinal Transfer Functions 

|U4|(deg/tec) _ / A 

w z(ib«) ^zCita) • “long 


Ngft’-Ngff/A long 


fcjq(deg/*c) _ 
^ZClbO 


2L36E-04 

(0) 
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(0.40) 

(0.61) 

2L33E-04 

(0) 

(L4E-06) 

(0.635 , 

0.563) 

L976E-04 

(0) 

(0.292) 

(0.644) 


L955E-04 

-(0)- 

(Z68E-06) 

(0.819 , 

0.508) 




















i gr(<kg) 
^(lb.) 


^long 


2.796E-05 

(0) 

(0.203) 

(0.370 , 

3.008) 

-L819E-05 

(0) 

(0.364) 

(2.255) 

(-4.452) 

" 2.130E-05 

(0.167) 

(0.351 , 

3.038) 


L470E-05 

(0) 

(0.261) 

(0.460 , 3.426) 



(L438E-05) 

(3.918E-02 , 

0.130) 

(0.652 , 

L382) 

(1.101E-05) 

(7.423E- 02 , 

0.147) 

(0.596 , 

L375) 

(3.949E-02 , 

0.118) 

(0.649 , 

L301) 


(L878E-05) 

(7.190E-02 , 

ai38) 

(0.588 , 

L279) 


Lateral Tra nsfer Functions 


KSS-Ngy/A, 


lat * 


N f (dC8) - N* (deg) / A 

^ZOta) *Z(lb«) ' “l 


lat > 


N0(<fcS) 

*^Z(lbi) 


^Z(Ib») 


-L58E-03 

(-.0805) 

(.927) 

-L59E-03 

(-.0922) 

(.904) 

-L43E-03 

(-.0723) 

(.981) 

-L44E-03 

(-.0879) 

(.940) 


3.19E-04 

(.468 , 3.65) 

2.15E-04 

(.611 , 4.17) 

2.89E-04 

(.447 , 3.96) 

2.04E-04 

(.593 , 4.33) 


(.0001) 

(L01) 

(.116 , L05) 

(.006) 

(L05) 

(.067 , 0.93) 

(.0028) 

(1.06) 

(.114 , L08) 

(.0065) 

0*09) 

(.060, .944) 


The engine transfer function for all configurations is given in short form 
notation by: 


C 


*(•»> 

*.<*> 


275 


(0.55X5) 



ation Storage Table for Quadruples 

*Pxxxx.4U : Quadruple with four engine inputs 
*Pxxxx.lU : Quadruple with one total engine input 
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